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Abstract 

The quantum periodic XXZ chain with alternating spins is studied. The properties of the 
related R-matrix and Hamiltonians are discussed. A compact expression for the ground state 
energy is obtained. The corresponding conformal anomaly is found via the finite-size computations 
and also by means of the Bethe ansatz method. In the presence of an external magnetic field, the 
magnetic susceptibility is derived. The results are also generalized to the case of a chain containing 
I different spins. 

CD . 

^ ■ 1 Introduction 

> : 

(1+1) -dimensional integrable - exactly solvable - field theories and one dimensional integrable lattice 
models have been extensively studied during the last two decades. They found numerous applications 
in a variety of areas in statistical physics ^ 12], condensed matter physics and high energy 
physics 0]. Apart from that, integrable models have their own mathematical interest in the context 
of quantum groups jSJ E] • 

The most known examples of quantum integrable models are the sine-Gordon model for which 
physical quantities such as spectrum and exact scattering matrices have been determined (Jj, and 
the Heisenberg model solved originally by Bethe [Sj. For the Heisenberg model as well, quantities of 
physical interest, such as spectrum, scattering matrices, and correlation functions have been explicitly 
derived 011131 • 

An integrable generalization of the Heisenberg model to higher spins and its g-deformation (the 
XXZ model) have been extensively studied in the literature (see e.g., ^2~GB1)- in this paper we 
will investigate the integrable alternating XXZ spin chain, i.e., an inhomogeneous chain with spin 
Si at odd sites and spin S2 at even sites. The first study of such a model, with (Si, S2) = (| 5 1) 
was done in [T^]. Later, various aspects of alternating chains were considered also in [20|-[2!)] but 
mostly for specific cases like (|, 1), (%, S), and (S, S+^). Our aim is to elaborate the generic (Si, S2) 
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case in detail in the framework of the algebraic Bethe ansatz with a particular emphasis on the 
thermodynamic limit and the finite size correction analysis and then to generalize the results to the 
case of I spins. In the present paper we will consider only the case for periodic boundary conditions. 
Some related results on open chains were presented in |26| . 

The paper is organized as follows. In Section 2 we outline the basic ingredients of the algebraic 
Bethe ansatz approach for the model in question, in particular, we discuss the structure and properties 
of the related R-matrix and Hamiltonians (technical details are given in the Appendix A). In Section 
3 we consider the model at zero temperature. We show that its ground state energy can be expressed 
in terms of a specific q-deformation of the ^-function (related to the Barnes double gamma function). 
Here we also present some finite-size computations of the corresponding central charge. In Section 4 
we first derive the Bethe equations in the thermodynamic limit at finite temperature T and in the 
presence of a magnetic field H. Then we compute the resulting magnetic susceptibility in the T — > 
limit. We also derive the effective conformal anomaly of the system by means of the thermodynamic 
limit of the algebraic Bethe ansatz and remark that the underlying conformal theory should have a 
spinon basis related to the D n algebra. In Section 5 we generalize our results to the case of the XXZ 
chain containing I different spins. 

For compactness of notations, we will use in the text both the deformation parameter 7 and q = e 17 . 
In Section 2 we assume that q is either real or takes values on the unit circle, whereas the results of 
Sections 3, 4, and 5 are valid only for the latter case. 

2 Algebraic Bethe ansatz 
2.1 R-matrix 

The underlying algebra of the XXZ model is the quantum Lie algebra U q {sl2) with the following 
defining relations 

[S + ,S-] = [2S%, [S 3 ,S ± } = ±S ± , (1) 

where 

is the standard definition of q-number. For generic q algebra (J2) has the same structure of represen- 
tations as the undeformed algebra sli- In particular, irreducible highest weight representations Vs 
are parameterized by non-negative integer of half-integer spin S and are (2S+ l)-dimensional. In this 
case, the generators (^Q) regarded as elements of End(C 2S+1 ) admit the following matrix realization 

(S + ) m ,n = {S~) n ,m = \J [m\ q [2S + 1 - m] q 5 m+ l,„ , (S 3 ) m>n = (S + 1 - 17l) 5 m , n , (3) 

where m, n = 1, 2S'-|-1. 

An L-operator associated with the XXZ model can be chosen as follows 




The defining relations Q are then equivalent to the following matrix exchange relation 

4? (A) Lf 3 (A + /x) Lf 3 ( M ) = L S M Lf 3 (A + M ) fl«(A) , (5) 
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where R^^(X) = La (A + |) (up to a normalization, see below). The lower indices in © specify 
components of the tensor product Vi®Vi (&Vs where the R-matrix and the L-operators act. Actually, 

2 2 

the R-matrix in © is just a particular representation of a more general object, R(X) G [[/^(s^)]® 2 ) 
that satisfies the Yang-Baxter equation 

flia(A) i?ia(A + #2 3 (m) = # 23 (^) i?i 3 (A + fx) flia(A) . (6) 

For a homogeneous spin chain or a related lattice model, it suffices to consider the "diagonal" 
solution R SS (X) of ©• However, to deal with an inhomogeneous chain, it is unavoidable to use the 
general solution R SlS2 (X). We provide the necessary details below. 

Define the "total spin" operator J on Vs 1 <8> Vs 2 as 

where (55 = |<Si — S^l, and Vj is the projector onto Vj in the Clebsch-Gordan decomposition 
of Vs 1 <g> Vs 2 . If f(x) is a function non-singular on sped, then eq. (J7J) allows us to define /(J) = 
Efijf 2 With these notations, a solution # 5l& (A) to the L r (? (s/ 2 )-Yang-Baxter equation ® 

consistent with the L-operator © can be written as follows (see Appendix A) 

R SlS2 (X) = (-iy- s i- s 2 ^(5i(Si+l)+S 2 (5 2 +l)-J(J+l)) x 

iyas + i + iA) iyj + i-iA) SlS ■ t 
x r 9 (5S + i-iA) r 9 (j + i + iA) ^ j ' e " ± - (8) 

Here r 9 (x) is a q-deformation of the gamma function such that 1 

T q (x + l) = [x] q T q (x), (9) 
and R± are the constant universal R-matrices [51 147] . 

~ g 4("-" 2 ) 

The g = 1 limit of (jSJ) gives the well-known s/2 _ solution constructed in |27j . 
As we discuss in Appendix A, the following relation holds 

(_1)J"25 q e(2S(S+l)-]Q+l)) ffiSSyl = p ^ 

where P is the permutation in Vs®Vg. Identity ensures an important property of the solution (|8|): 

i? 55 (0)=P. (12) 

Also, one can show that (jSJ) enjoys the following symmetries 

R 1/q (X) = (R.i-X))- 1 = (^(A))* , (13) 
R 1/q (0) = (RM)- 1 = (R^O)) 1 = R q (0) , (14) 

where t stands for the matrix transposition consistent with (J3J). 

It should be stressed here that eq. © alone does not define the q-gamma function uniquely. However, we need to 
deal only with expressions of the form T q (x + n)/V q {x), n G N that are well defined. 



R f 



Errir (^-^^Y^ 3 . ao) 
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In order to obtain R SlS2 (X) in terms of the algebra generators, which is needed, e.g., for construct- 
ing related Hamiltonians, one can substitute into (jSJ the following explicit form of the projectors (see, 
e.g., ES|) 



Si+s 2 

n 

l=5S 



x s lS2 _ [l]q[l + 1]q 

[ 3 -l] q [j + l + l] q ' 



(15) 



where 



X s ^ = (q s3 ® 1) (5+ <g> S~ + S~ <g> S+) (1 <g> q~ s3 ) + 



S 3 N 



cos 7 
2 sin 2 7 



'1 <g> 1 + g 253 ® g" 253 



,. .— 1 (1 <g> q 

2 sin 71^ ^ 



-2S 3 



cos7(2,Si + 1) + (g 2S ° 1) cos 7 (25 2 + 1' 



(16) 



Formulae (jHJ), (fT31) . and (fTB|) allow us, in particular, to verify that Rz s (X) does reproduce the 
L-operator Let us however describe another way of doing that. To this end we will construct the 
"Baxterized" form of -R(A) for higher spins. As we show in Appendix A, the projectors admit the 
following form 



S1+S2 g 2l(l+l) _ 9 2Si(5i+1)+2S 2 (5 2 +1) jR -1 jR _ 
Vj = II g 2i(Z+l) _ q 2j(j+l) 



(17) 



i=<5S 



This formula along with (jHJ) allows us to find an expression for R(\) in terms of -R+ and R- : 

(q-q- 1 )- 1 



R^ S (X) 



R 1S (X) 



qh~ iX (El S )- 1 -^ x -h(Ri S )- 1 



2 1 " K \<i 
(q-q- 1 )- 2 



(18) 



[S + iA] g [S + 1 + iA] 



q 



2iA-l 



q~ z w ) (i^ 5 )" 1 + <? 4 tu (R^y'R^iR 



r>lS( olS\-l 



+ 



1-2iA 



+ (l-g 2 ) ro - (? 2 ^ 1 )(i?f)- ] 



[25 + 1], 
[4S + 2] 9 ! 



(19) 



etc. From ()18[) one easily verifies that [S + 1 + iX] q R2 S (X — i) coincides with L S (X) given by (JIJ. 

Let us remark that, twisting we can obtain other L-operators related to the XXZ chain. In 
particular, the following one is often used in the literature: 



L(A) = e 



7 (A-i)S 3 



L(A)e 



7(4"A)5 3 



[5 3 - iA] s 



s- 



S + -[S S +iX]g 

From the Yang-Baxter equation (jHJ) and the symmetry of the solution Q 

[R ab (X),S 3 a + S 3 b ]=0 1 
it follows that the universal R-matrix corresponding to (f2*0|) is related to © as 

It has the following symmetries (see Appendix A) 



(20) 



(21) 



(22) 



R(-X) = (R(X))- 1 , R 1/q (-X) = R q (X) , {R(X)f = R(X) 



(23) 
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2.2 Hamiltonian 

In order to construct a Hamiltonian for a closed alternating spin chain we apply a generalization of 
the method developed in ;29 and consider two transfer matrices: 2 

r«(A) =tr a TW(A) = trjR s a f(\) R s a f^(\) . . . R S J>(\) R s a f(\)) , 1 = 1,2, (24) 



where the auxiliary space (denoted by the subscript a) is Vsi and V£ 2 , respectively. Since the mon- 
odromy matrices satisfy the Yang-Baxter equation 

R S /i (A - n) if (A) if Oi) = if (//) T« (A) JZ* S ' (A - M ) (25) 

the corresponding transfer matrices commute 

[r«(A),r«)( At )] = 0. (26) 

Therefore, the first and higher order logarithmic derivatives of rW (A) at A = are quantum integrals 
of motion. Moreover, thanks to the property IJ12JI . these integrals are local, according to Theorem 3 
in [SOI. 

Using the property (|T2*|) (along with tr a P a t = 1), we find that 

AT-3 AT-l JV-2 N 

r«(0) = ( n P M+2 ) ( J] /$&(())), r( 2 )(0) = ( ]J F fc , fe+2 ) ( ]J ^(0)), (27) 

fc— odd fc— odd fc— even fc— even 

where iV + 1 = 1 (the periodic boundary conditions). Noticing that for both r(0)'s in (|27|) all R(0)'s 
commute among each over, it is not difficult to obtain the following expressions for Hamiltonians 

N-i . N 

= ia A lnr«(A) H [ n\ W W = i^lnr^A) £ (28) 



A=0 



A=0 

n— odd n— even 



(i) 

with Hn being local Hamiltonians involving three nearest sites: 



Hj? = (<5-i(0)J dxR^iX) a=o + 2 ^i(O) , (29) 

where 1 = i+1 (mod 2), and 



2 



A=0 



is the local Hamiltonian of the homogeneous spin S chain. These formulae provide a particular case 
of local Hamiltonians for inhomogeneous chains derived in [3U] (see also ^1 ) ■ Obviously, they recover 
the homogeneous case as well since for S\ = S2 we have R^n+ii®) = ^Vn+ij an d then H$ becomes 
a sum of two terms of the type (|3U|). 

Let us notice that, computing in the same way the local Hamiltonians corresponding to the 
R-matrix (1221) and taking the property (|21|) into account, we obtain 

H^=H^ + l ^(Sl-Sl +2 ). (31) 



2 If we introduce proper inhomogeneities of the spectral parameter, we will obtain massless relativistic dispersion 



relations for the particle-like excitations of the model |2o|. 



5 



Hence for the total Hamiltonians we have 7i« = H^. 

Let ^ q {x) denote logarithmic derivative of the g-gamma function used in Q, so that we have, in 
particular, *$> q (x + 1) = ^ q (x) + 7 cot 72;. Then jSJ) allows us to rewrite the local Hamiltonian in the 
following, more explicit, form 

F« = (R e )n,n+1 q^+^, n+ x+l) (_!)!«,«+! [* ? (J n) „ +1 + 1) - * q (5S + 1) 

+ ^(Jn,„+2 + 1) " (-l) 1 "^ 1 s -d»,n+i(Jn,n+l+l) ( J R £ - 1 ) n>n+1 , (32) 

where the quantum space is V^v at the sites n, n+2, and Vs ? at the site n+1. 

In order to exemplify (|29|) and (|32|) . let us construct explicitly for (Si, S2) = (|, 5). Consider 
the XXX case first. Here we need only the following projectors 

C^km = + \ ± \ ± ^ ' ' (33) 

where <r = (cr , a 2 , a 3 ) consists of the Pauli matrices. In this case (|3*2*)l yields 

+ {(<r n • S n+1 ), (5 n +i • ^n +2 )} + (\~ 5(5 + 1)) (<?„ • 0W2) + |) . (34) 

Here and below { , } stands for the anticommutator. Since {(<? • a), (a ■ b)} = 2(o ■ b), it is easy to see 
that, for 5 = 2> ec l- (EU reduces to the usual Heisenberg Hamiltonian. In a slightly different form 
(|34|) was given also in |2U| . 

In the XXZ case, it is easier to compute Hn\ taking into account the symmetry Q23JI . To write 
down the answer in a compact form, we introduce the Sklyanin generators |31j of Uq^sfo), 

&° = [\} q co S ( 7 S 3 ), S 1 = i(5+ + 5-), S 2 = l(5 + -5-), 6 3 = cos(^)[5 3 ] 9 . (35) 
They satisfy the following quadratic relations 

[S 3 , 5±] = ± {S°, 5 ± } , [6°, 5±] = T (tan ^) 2 {6 3 , 5 ± } . (36) 

Using these relations, we obtain from (|29|) - (|31|) the following (i, 5)-Hamiltonian (see also 19 for a 
(2, l)-Hamiltonian) 

^ = Alq I 112 ( 2 { 6 n+l> ( g n ' + (®n + l ' <% +2 )}+ COS^S*) 

+ {(<? n • 6 n+1 ), (6n+i • al +2 )} + ((6° +1 ) 2 - 6 n+1 • S n+1 ) (a n • 3% +2 ) 

+ j cos 7 + [6° +1 , [{a n ■ 5 n+ i), (5„+i • ct^ +2 )]]) , (37) 



where 6 = (S 1 ,© 2 ,© 3 ) and a 7 = (a 1 , a 2 , a 3 cos 7). For 5 = | we have 6° = |, & = \B and then 
(|37|) reduces to the usual XXZ-deformation of the Heisenberg Hamiltonian. 

Consider the *-structure on Uqisl-i) corresponding to the compact real form U q (su(2)), i.e., an 
anti-automorphism such that 

(S^)* = S T , (5 3 )* = 5 3 (38) 
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and that extends onto a tensor product as (£ <g> £)* = ^* (g> £*. If g is real or |g| = 1 we have (see 
Appendix A) the following "unitarity" relations 

(R q (\))* = (i?q(A)) _1 , (R q (\)Y = (^(A))" 1 • (39) 

Applying the latter relation along with (|14[) to the counterpart of (|29j) for .ffj^ and taking into account 
that h* nm = h nm (see, e.g. [2H1), we infer that 

{H®Y = fl#> . (40) 

Therefore, the total Hamiltonian of an alternating spin chain Ti® is hermitian if 7 E R or € R. 

2.3 Momentum operator 

Unlike the homogeneous case, the transfer matrix r®(0) is not the shift operator and hence its 
logarithm cannot be interpreted as the momentum operator. Observe however that the products of 
permutations entering (|27j). 

N-3 N-2 

= n f ^+2 , = yi Ffc)fc+2 , (41) 

A;— odd k— even 

are shift operators of odd/even lattice sites by two steps, i.e., 

0-0 £ n = £ n+2 tf-M if n = i mod 2 . (42) 

This motivates to consider quantum integrals of motion generated by r 1,2 (A) = (A)r( 2 ) (A). Due 
to (|26[). the corresponding Hamiltonian is given by 



, J 
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J 



W = 49 A lnr( 1 ' 2 )(A) = ^ (W« + W< 2 )) (43) 



A=0 2 



with H® defined in and .7 > being a coupling constant. 

In order to compute lnr 1,2 (0) we need a property of the R-matrix that follows from evaluating 
the Yang-Baxter equation © in Vs 1 O Vs 2 <g> Vs 1 at /U = — A. In view of (|12|1 this leads to the relation 

Rf^(X)Rlf^-X) =R s 2 fi(-\)R s 3 ^(\), (44) 

and we infer, in particular, that 

i?f 2 l52 (0) Rlf 1 (0) = 1 1 (45) 

if the normalization of the R-matrix was chosen appropriately (as, e.g., in |5J)). 
Now, using ip7|). (g2J), and (gSJ), we obtain 

T i> 2 (0) = £/« (46) 

which generates a shift by two lattice sites. Thus, V = A hiT 1,2 (0) can be regarded as the momentum 
operator (as was observed also in ^H]). Similarly to the homogeneous spin chain, the pair (Ti,V) 
defines in the gapless regime (\q\ = 1) at zero temperature a conformally invariant system. 
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r (0) 



2.4 Bethe equations 

Consider the alternating XXZ chain of even length JV. According to the general scheme of the algebraic 
Bethe ansatz |291 I3U| (see also [32] for a review) , eigenvectors and eigenvalues of the auxiliary transfer 
matrix 

(A) =tr a Tj°>(A) =tr a (4 2 (A)4 1 _ 1 (A)...L^(A)Lf 1 (A)) , (47) 
where the trace is taken over the auxiliary space Vi , are parameterized by solutions of the system of 

2 

equations 

2 — M 

where M is the level of a Bethe vector \I/(Ai, . . . , Am); it is related to the total z-component of the 
spin as 

N 

(J2 Si) • • • , A M ) = S z *(A X , • • • , A M ) , (49) 

n=l 

where = ((Si + S 2 )N/2 - M). The Yang-Baxter equation 

(A - /x) Ti°) (A) if (/i) = if (//) Ti°) (A) fljf (A - M) , j = 1, 2 (50) 

implies commutativity of the corresponding transfer matrices, 

[t(°>(A),t0>(m)]=O. (51) 

Therefore the eigenvectors of r^°'(A) are simultaneously eigenvectors for r^(A). Moreover, as follows 
from the results of jHU], eigenvalues of the momentum operator and the Hamiltonian Tt (in the 
presence of an external magnetic field H) are given by 

M 2 j M 2 

P(A X) . . . , Am) = £ X)p(Aal-Si) , E(X h Am) = --^ 5>'(A a |^) - iTS* , (52) 

a=l j'=l a=l j=l 

where we introduced the functions 

p(X\S) = — In — — — -r = arctan(tanh(7A) cotHS)) , (53) 

2i smh7(A + io) 

The corresponding values for the alternating XXX chain are obtained by the 7 — ► limit of (|52 |) -()54 |) . 

3 Ground state energy and finite size effects at T = 
3.1 Ground state energy 

Let ir/j be integer. In the thermodynamic limit, when M, N — > 00, the string hypothesis holds |33| . 
namely solutions of the Bethe equations (|48|) form string-like clusters, i.e., groups of the type 

Xt k) = X n a + l -{n + l-2k), fc = l,2,...,n, (55) 



A (°) = A°+i— 
' 27 ' 
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where the centers of strings A™ and A^ are real, and A^ ) is the negative parity string. The number n 
is called the length of a string. 

At the zero temperature the system is in its ground state which, as we will see below, is a Dirac sea 
filled by strings of lengths 2S\ and 2£ 2 (negative parity strings are not present). Taking logarithm of 
(|48|) and carrying out summation along the strings in the standard way, we obtain a set of equations 
on the centers of these strings: 



N 



Ep(Ag ) N+ E P^HJ =^ ) + EE^,25 J (A«-A^ ) ), (56) 

m=i m=«5S+i J=l/3=1 



where (55 = |Si — S2I and Aq^ = A^ 5i , i = 1,2. The numbers are distinct integer or half-integer 



numbers such that \Qa^\ < y — \, and 



^(n+m) — 1 

*n,m(A)= E (p(A|0+P(A|i + l)) (57) 

!=2 I"— rn\ 



(for n = m the term with / = is omitted). 

When N,M — > 00, ()56j) become integral equations: 

2Si — I S1+.S2 — 5 2 



£p'(A|m)+ E P / (A|m) = 2vrp«(A) + 2^(cI>' 2Ssi25j * /3 W)(A), (58) 



m=| m=5S+\ J =1 



where the convolution is defined as {u * v)(X) = dfiu(X — fi)v(n), and <fr' n m (A) is the derivative 
of (|5T|) . It is straightforward to check that 



d\e lkX p(\\m) = 7T if m < '- ■. (59) 

27 



sinh(g - mfc) ^ 

! — r if m < — 

sinh ^ 7 



Using this and taking Fourier transform of (|58|) . it is not difficult to see that the only solution of 
these equations is given by 

P W = P {2) = 1 if 2 7 max(5i , 5 2 ) < vr . (60) 

4 cosh 7rA 

As follows from (|52*|) , (|54")l , and (|6T7)) , computation of the ground state energy amounts to computing 
a sum of integrals of the following type 

T( x r P'(Mm) f°°„ 1 sinh(^-mfc) 

I 7 (m)=/ dA 7-^r= dk r -J- £ - , 61 

coshvrA J cosh § smh f£ 

where, in the last equality, we took ()59|) into account. Expanding the functions in the denominator 
in series, we obtain 

n>l 



m 3 irn\ / m 1 irn 



2 4 27/ V 2 4 2 7 



(62) 
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Let us now show that 7 7 (m) admits, for real 7, a compact form in terms of a certain deformation 
of the ^-function. To this end we define the following special function 



r (x) - (q _ g -i r * r 2 (bx\b,b i) 



(63) 



where q = e 17rfe2 , 0<6<1 and r 2 (a;|a;i, 0*2) is the double gamma function introduced by Barnes [3~I] . 
Properties of r2(x|wi, W2) imply the following relations 



F q (x + 1) = [x) q T q {x) , (5 - q- l ) l,b2 T q {x + l/b 2 ) = 2 sin(vrx) F q (x) , 



(64) 



first of which allows us to regard T q (x) as a ^-deformation of the gamma function. T q (x) is closely 
related to the noncompact quantum dilogarithm that has been actively studied recently [33] in the 
context of integrable models. In particular, it can be shown that 



hxT q (x) = const -xln(g-g 1 ) + ^(b + b 1 - bx)bx - J y (1 _ e te)(i _ e «/6) " 



M+iO 



(65) 



This expression allows us to establish connection between the integral in H61|) and the g-deformation 
of the 'I' -function defined via ^ q (x) = d x lnr g (x). Namely, we find that 



m 3 



J 7 (m) = ^ q 2 — + - - 2 — + - , where q = e 



m 1 



=17 



(66) 



Now, from (|52|) . 1)60(1 . 1)61(1. and 1)66(1 . we obtain the ground state energy per lattice site 



J 
4 



e(Si,S 2 ) 



J 



2 2Si-i 



Sl+S2 — -R 



»=1 m=i m=<5S+| 
J " 



(67) 



For a homogeneous chain, (|67j) reduces to 

Jr 



e(5) 



2 L 9 



M 5+ 2 



(68) 



Replacing q-deformed \& -functions with ordinary ^-functions in (|67 j) -()68 |) . one recovers the energy 
densities in the XXX case. In particular, (|68j) turns in this limit into the well-known formula |12j . 

It is worth noticing that the first relation in (|64[) implies that ^qi{x + 1) — ^ q i{x) = 27cot27x. 
Therefore, if S\ and S 2 or, respectively, S are integer, we can rewrite 1)67(1 and ((68(1 in terms of 
trigonometric functions: 



S2-I 



5(5i+Sa)-l 



r(.S',.5 2 ) --^(^cotfrpfc + l)] + V eo([;(2A- + .1 )] +2 V cot 7 ! 2/,- + I. )] ) . (09) 



k=0 k=0 
S-l 

e{S) = - 7 J £ cot[ 7 (2A: + 1)] . 

fc=0 



fc=±<5S 



(70) 
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3.2 Finite size corrections 



It is well known that statistical systems at criticality are expected to exhibit conformal invariance |3b| . 
Therefore, the critical behavior of such systems should be described by a certain conformal field theory. 
One can identify the corresponding effective central charge c by investigating the finite size effects 
of the ground state of a critical system [37j. For a spin chain of a finite length N, the ground state 
energy at zero temperature depends on N as follows: 

e N = e 00 -^ + o(N-"), (71) 

where v s is the speed of sound. With the energy and momentum normalization as in (|52|) we have 
v s = ttJ/2 provided that the condition in (|rJU)) is satisfied ^3E3I1I]- I n this section we choose J = 1, 
so that v s = 7r/2. 

We will apply formula (|71|) to the alternating XXX spin chain in order to find the corresponding 
central charge numerically. To this end we start with solving numerically eqs. Q56JI in the case 7 = 0. 
The solution of these equations is unique since they are extremum conditions for a convex function 
(Yang's functional, see |2j for the analogous consideration in the homogeneous case). 

It is interesting to notice that, if we assume that roots of the Bethe equations form exact strings 
even for finite N and that the centers of these strings are given by (|56|) , then formula Q71|) yields the 
same value c = 2 + 0(A r ~ 1 ) for all alternating chains (see Table 1) thus predicting the conformal 
anomaly c = 2 independent of the spins (provided that S\ ^ S^)- The interpretation of this phe- 
nomenon is that the field theories arising as H = 0, T — > limit and T = 0, H — > limit of a spin 
chain are, in general, different (see the discussion in |22| ) and a finite-size computation based on the 
string hypothesis corresponds to the latter limit. The same phenomenon occurs for a homogeneous 
chain; in this case an analytic or numeric finite-size computation using the string hypothesis gives 
c = 1 independent of spin S, whereas those not using the string hypothesis give c = 3S/(S + 1). 
Comparing this fact with our result described above, we may conjecture that a finite-size computa- 
tion based on the string hypothesis performed for a chain containing equal number of I different spins 
will give c = I (which is the smallest possible value of c for such a chain, see Section 5). 



Si, S2 




A S2 


c 




c 




±0.129530 
±0.538739 


±0.131152 
±0.557922 


2.08 


±0.543466, ±0.129918, 
±0.547454 ±i 0.504109, 
±0.130894 ±i 0.500000 


2.06 


1,1 


±0.130724 
±0.554788 


±0.131252 
±0.561180 


2.10 


±0.549334 ±i 0.552225, 
±0.130456 ±i 0.516086, 
±0.562630, ±0.561099 ±i 1.007705, 
±0.131684, ±0.131684 ±i 1.000000 


2.59 


1 3 
2' 2 


±0.130351 
±0.535032 


±0.131852 
±0.552153 


2.08 


±0.449421, ±0.094976, 
±0.668491, ±0.544016 ±i 1.055667, 
±0.168193, ±0.131555 ±i 1.017436 


2.60 



Table 1. Alternating XXX chains with N = 16 sites. A^ are the real roots of eq. I|56|l for 7^0; the estimate 
c of the central charge is computed by (|71J) under the assumption that Af i are centers of exact strings. A™ are 
the complex roots of eq. (|48|l for 7 — > corresponding to the genuine ground state; c is the corresponding 
finite-size estimate computed by (I7T1) . 

Although, as we discussed above, the string solutions constructed on the base of (|5tj)) are not 
exact, they, as seen from Table 1, provide good initial data for iterative numerical computations. 
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Using them as a first approximation and applying a Newton- type method, we solved numerically the 
XXX version of eqs. (|48[) for several alternating chains and computed the corresponding ground state 
energy. Then we used formula (|71|) and eq. (|67[) in order to find the finite-size estimate of the central 
charge for these chains. Our results, presented in Table 2, are described by the following expression 

_ 35i 3{S 2 - Si) , 
C " S 1 + l + S 2 -Sx + l- {U) 

This formula was conjectured in [2Q' a as an extrapolation of the results for the (o,S) XXX chain. 
Since it is known that for the homogeneous XXZ chain the central charge does not depend on the 
anisotropy (if 275 < it) and coincides with that of the XXX chain, one can expect that (|72|) applies 
to the alternating XXZ chain as well (provided that the condition in Q6U|) is satisfied). This is indeed 
so, as we will demonstrate below. 





iV = 8 


N = 16 


N = 24 


N = 00 


Si, 5*2 


|e| 


c 


|e| 


c 


|e| 


c 


|e| 


c 


I>1 


0.666301 


2.21 


0.644511 


2.06 


0.640791 


2.03 


0.637888 


2 


Li 


0.917421 


2.79 


0.889926 


2.59 


0.885257 


2.55 


0.881620 


2.5 


2 2 

2' Z 


1.086631 


3.14 


1.055557 


2.91 


1.050305 


2.86 


1.046222 


2.8 


2 & 
z ' 2 


1.213910 


3.39 


1.180345 


3.12 


1.174695 


3.07 


1.170311 


3 


1 3 

2' 2 


0.661489 


2.84 


0.633361 


2.60 


0.628649 


2.56 


0.625 


2.5 


i 2 

2' Z 


0.666989 


3.24 


0.634793 


2.94 


* 




0.625352 


2.8 


1 5 
2' 2 


0.675916 


3.52 


0.640833 


3.17 






0.630647 


3 



Table 2. Alternating XXX chains with N sites (the coupling constant is J = 1). e is the energy per site for 
the numerically found ground state solution of eq. I|48[l for 7^0; c is the corresponding finite-size estimate of 
the central charge obtained by eq. J7TJ. The N — 00 values are predicted by eqs. Q67f) and l|72|) . The asterisks 
indicate lack of stable numerical results. 



4 Thermodynamics 

In this section we will investigate the thermodynamics of the alternating XXZ chain in the case where 
the anisotropy parameter v = — > 3 is an integer number and v > 2Si. Our analysis follows, with 
proper modifications, the standard Bethe ansatz method [9l ll"4"| UH l H9 l 132]. 

4.1 Bethe equations at T > 

Denote i9j = 2Si] without loss of generality we will always assume that $2 > Following |19| I25j. 
it is convenient to introduce the functions 

^_ sinh 7 (A + f) iTr cosh7(A + f) 

e " (A;i/) = smh 7 (A-f)' 5 " (A;i/) = e " (A± 2^ ) = cosh7(A-f)- ^ 

Then the Bethe equations ()48|) acquire the form 

M 

e^A^e^A^ 2 = - J] e 2 (A Q - A^) . (74) 

13=1 
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The spin, energy and momentum of a state are given in terms of the Bethe ansatz roots X a by eqs. (|49|) 
and (|52[). In this section we set J = 1/tc. 

In the thermodynamic limit N, M — > co the string hypothesis holds , and the Bethe equations 
become 

2 u—1 M m M 

\{X n ^{\l) N l 2 = {-) n \[J[E nm {\l-W[[G nl {\l-\%, (75) 

3=1 m=lf3=l 0=1 

2 v-1 M m M 

n^-^)^ 2 = - n n - ^ n ^ - ^ . ^ 

i=l m=l/3=l /3=1 

where the notation for the centers of strings are as in (|55|). n takes values 1, 2, . . . , (y — 1), and 

-Xnm(A) = e| n _ m+1 |(A) e| n _ m+3 |(A) . . . e( n+m _ 3 )(A) e( n+m _ 1 )(A) , 

E nm (X) = e.\n-m\ W ef n - m+ 2\ W ■ ■ ■ W e {n+m) W , 

G n m{\) = 9\n-m\W 9\ n - m +2\( X ) ■ ■ ■ 9{n+m-2)W 9{n+m)W ■ (77) 

The energy (|52|1 of a given state takes the form 

v— l M n n Mo 



where 



n=l a=l a=l 



Z nm (X) = -^-—lnX nm (X) , b n (\) = -^-—\ng n {\) . (79) 

Z7T UA Z1T UA 

It is useful to notice that the Fourier transforms of these quantities are 



sinhj^ — max(n, m))^J sinh^(min(n, m))^ 
smh(— )smh(|) 



b n (u;u) = 0<n<u, = V , v<n<Sv. (81) 



sinh (m sinhf(n-2i/) 2 

- . - )?' , 0<n<u, = Vtt 

sinh(^) sinh( 

In the N — > oo limit we obtain from ()78[) an expression for the energy density: 

v-l 



E 

N " 2 



1 ' ^ poo 1 /*oo 

-J^/ d\{Z n ^{\) + Z n ^{X)) Pn {\) - - d\(b^(X) + b^(X)) Po (X) 

f 1 /*oo /*oo 1 

+#(Vn/ d\p n (\)+ dApo(A)--(^i+^)), (82) 

where p n is the density of the n-strings (pseudo-particles) and po is the density of the negative parity 
string. The energy density of the ground state (T = 0, H = 0), which consists of two filled Dirac seas 
(strings of length i?i and $2) is given by 

e o = ¥ = -lE dAZ* 4 *,(A)fl(A), (83) 
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where (cf. ifSO))) 

z cosn(7rA) 2 cosn ^ 

Besides the densities of pseudo-particles we need to introduce the densities of holes p n (for n-strings) 
and po (for the negative parity string). Then, taking the logarithm of l|75 j) -([76 j) and differentiating, 
one finds 

1 u ~ 1 
Pn(A) = -{Z n , &1 (X) + Z n $ 2 {\)) - ^2 A nm * p m (\) - B n * po(X) , 

m=l 

1 V ~ l 

-(Po(A) + poW) = 2 + ^ 2 ) - 5m * ^ m ( A ) ~ ° 2 * ' ( 85 ) 

m=l 

where 

Aim (A) = (5 nm (5(A) + ^-—lnE nm (X) , £?n (A) = 77- 4r In G\ n ( A) , a n (A) = lne n (A) , 
Z7T aA Z7T aA Z7T aA 

2 cosh( — ) sinh(— ) 

A nm {u) = 2cosh(-) Z nm (uj) , B n (uS) = . \ fVLJ . — h ^.f-i • (86) 

In order to find p n , po in terms of p n , p~o, it is convenient to consider the convolution of the first of 
the density equations (|85[) with the inverse of A nm which is given by 

A nm = 5nm ~ s( U )(^n,m+l + K,m-l) ■ (87) 

Taking into account the following identities (summation over m is implied) 

A nm* Z mkW = S{\)5 nk , * 771 = § 5 n ^_l , ^ * B m (X) = s(X) 5 n>i ,-2 , (88) 

a 2 (A) + s*^_ 2 (A) = 0, b n (X)+s{X)5 n ^ 1 = -s*Z n ^ 2 (X) = s*B n (X), (89) 
we obtain the following expressions 

1 " _1 
Pn(A) = gg(A)(<Wi + <?n0 3 ) ~ / , 4l * Pm(A) + S n>v - 2 S * /Pp(A) , 

m=l 

Po(A) = -po(A)+s*p v _ 2 (A) + -s(A)(<5 1?1)I ,_i + (^ 2)1/ _ 1 ). (90) 
According to the standard statistical principles, the entropy density of the system is given by 

s = I = £ r ^A(p n (A) ln(l + M^T) + ~ Pn {X) ln(l + ^)) . (91) 

r^nJ-00 V Pn(A) p„(A) J 



N 

n=0 ' 



Denote r/ n (A) = g4^. The free energy of the system is F = E — TS. The equilibrium condition, 
5F = 0, along with relations (|82jl. (|8ft|) . and (|9*T|) yields the thermodynamic Bethe ansatz equations 

Tln(l + r ?n (A)) = nH-~(Z n # 1 (X) + Z n # 2 (X)) 



+T J2 A nm * In (l + VmW) ~ TB n * hi (l + % '(A)), 



m=l 



r M * + -j, A ,l ) = ir-i(MA) + MA)) 



i + %W 2 



+T B m * In (l + ^(A)) - Ta 2 * In (l + r^\\)) . (92) 



m=l 
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These equations can be rewritten with the help of (|87|) - (|89l) as follows 

uH 

e n (A) = s(A)*Tln[(l + r ?n+1 (A))(l + r /n _ 1 (A))] + — 

+ S n> „- 2 S (A) * Tin (l + 77o X (A)) - ^s(A) + <5 n ^ 2 ) , (93) 

uH 1 
eo(A) = —-s(A)*Tln(l + ^_ 2 (A)) + -s(A) (^^1 + ^,^1), 

where e n (A) = Tlnr/ n (A) are the so-called pseudo-energies (the terms with ln(l + 770) and ln(l + r/ u ) 
are omitted in the equations for t\ and e v -i, respectively). Observe that the equations for eo and 
e v -% in (j93|) lead to the relation 

eo(A) = */#-€„_! (A). (94) 

This implies, in particular, that the case #2 = v—\ (recall that i?i < #2) is special because then an 
extra term appears in the equation for eo- Indeed, if $2 = ^— 1, the -^-strings get close to the odd 
parity strings, namely Q(A^°' ) — A^ 2 ' 1 ^) = 1 (see (|55|1). As seen from the r.h.s. of ()48[). this affects the 
structure of the roots of the Bethe equations. 

Using (|85j). (|92|). and (|89|) . we obtain the following expression for the free energy density at the 
equilibrium 

v-l 



f = ^ = -oE rfAln(l + r ? - 1 (A))(Z^ 1 (A) + Z ni?2 (A)) 

n=1 J-oo 

+ I jT dX ln(l + % 1 (A))(6^ 1 (A) + b» 2 (A)) - ( ^ + ^ )H 



-00 
2 



= eo-«Z)l/ dAs(A)ln(l + Wi(A))+^-i / dAs(A) ln(l + ^(A)) } . (95) 
4.2 The T = 0, if < 1 limit 

At zero temperature, the ground state in the anti ferromagnetic regime consists of two filled Dirac 
seas (strings of length #1 and $2 > In this case we can rewrite the equations (|9*2*)l for the strings 
#1 and $2 only in the following approximate form 



*(A) = ~\ E ^ ( A ) + #iH + J2 M*i * Tln ( X + <( A )) . ( 96 ) 
i=i j=i 

^■nm(A) = ^4 n m(A) — 5 nm 5(X) . 



In the T = limit the following quantities are defined for i = 1, 2 



Tln(l + r^ 1 ) - ±ef , e± = ^(e* ± . (97) 



Then equation (|96f) can be rewritten as follows 



2 2 2 

* ej(x) = -^J2 ( A ) + * H + E ^ * 4( A ) • ( 98 ) 



2 

i=i i=i 3=1 
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By solving the above system with the help of (|88|). we obtain 

e i (\) = -ls(\) + h i2 tH + M ij *eU\), £ = — V — , (99) 



where the kernel M has the following Fourier image 

hi{u) 

h(uj) h 2 {co) 



M(u) = ( W Hf\ ) , (100) 



with 



sinh(&?-l)f sinh(??i-l)f 
^ ~ 2 cosh | sinh + 2 cosh § sinh(t?i)f ' 



2 oiiihul> 2 ^won 2 oiij.ii ^ u i ; 2 

sinh(<5tf-l)f sinh(i/-tf 2 - l)f ?, ^_ sinh' 2 ( 

2 ^ ~ 2 cosh f sinh <5t?f + 2 cosh § sinh(zy - i? 2 )f ' ~ 2 cosh f sinh ' (H)1) 



and 5$ = $2 — Notice that, as in jHH], the quantities hi, h are related to the LL (RR) scattering 
between left (right) movers, and the LR (left-right movers) massless scattering (see e.g. [39 ). 

Eqs. (|99|) imply that £j(A) are monotonically increasing functions of |A| with the maximal values 
given by 

£l (oo) = ((1 - M(O))- 1 Q) . = ^ . (102) 

Therefore, for H small but non- vanishing, the pseudo-energies £j(A) are negative only on finite inter- 
vals |A| < aj. Having (|99() in mind, we introduce the zeroes of the pseudo-energies as follows 

ai = —(ln£H + lnki + 0(l/]nH)), (103) 

where ki (which do not depend on H) are to be determined later. Next we define the functions 

Si{\) = {£H)- X ei(\ + oti) for A > (104) 

and (A) = for A < 0. Making the shift A -> A + on in (JHU) and using that s(A + 04) £H e~ nX 
(since on — > 00 as i? — > 0), we obtain a Wiener-Hopf type system 

fci Z" 00 f 00 

5i(A) = — -e' wX + dX'hxiX- X')Si(X')+ / dA'/i(A - A' + t*i - a 2 )5 2 (A'), 

S 2 (A) = ^~4 e ' 7rX+ d\'h(X-X' + a 2 -a 1 )S 1 (X')+ dX'h 2 (X - X')S 2 (X'). (105) 
^ ^ jo io 

The above system can be solved via the standard Wiener-Hopf methods as follows. The Fourier 
image of the matrix kernel K of the system (|105j) is 

ire. ,\ — f e-Mai-oa)^) \ 

KM " I e- ; -(— h 2 (u) ) ■ (1 ° 6) 
Its resolvent admits the following factorization jlU] 

(l-k(u))- 1 = G+(u)G-(u), (107) 
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where G±(u>) are analytic for ±Qu> > 0, respectively, and also G\_{oj) = G_(— u) for 9u; > 0. The 
solution of the system ()105j) (see also (U) is 

2 UJ + lO CJ + 17T \1 J 

k = Gz\-m)G-(0)P\ . (109) 

Denote D = det(l — if (0)). From (|1U6|) — (|1U9|) we compute the following quantity 

tfSim) = -L( l)G + (0)G_(0)(J) = ^(1-^(0)) = ^. (110) 

Recall that the odd parity string contributes to the free energy (|95f) only if $j = i/— 1. Observe, 
however, that its contribution vanishes as T — > because it follows from (|102|) and (|94|1 that eo(A) is 
positive. Thus, as T — ► 0, the free energy of the system becomes 



/oo " poo 

dX s(X)ei(X) =e -^Hj2 dXs(X + a^(A) 
*=± - i=i Jo 

2 /■OO 

= e -( 2 H 2 J2 k i dXe-^SiiX) = e -£ 2 H 2 k t S(m). (Ill) 
i=i ^ 

Substituting here (|110j) . we obtain the following expression for the magnetic field dependence of the 
free energy, 

(112) 

2ir v — vi 

The magnetic susceptibility is given by x = ~ ^nf- The va l ue °f u sX does not depend on the 
normalization of the Hamiltonian (v s is the speed of sound, which in this section is |). Thus, we infer 
from lfTT3|) that 

v S X=o 5-- ( 113 

2tt v — vi 



Notice that the contribution of the external magnetic field to the free energy (|112|) does not depend 
on 6$, that is why our result coincides with the magnetization of the homogeneous XXZ chain found 
in |14j (where v s = Jtt/2). Also, choosing #2 = 2 in (|113|) . we recover the results obtained for the 
alternating (|,1) chain both in the XXX [23] and XXZ [2U cases. 

4.3 The H = 0, T < 1 limit 

Here we will investigate the low temperature thermodynamics of the model in order to find the 
corresponding conformal anomaly. In order to evaluate the entropy of the system at small but 
nonzero temperature we use the following approximations which hold for A> 1 (see, e.g., |141 116j ) 

p n (X) = I(_l)'»o /n(A) * (A ) , ~ Pn{X ) = i(-l) a «°(l - / n (A)) i n(A) , (H4) 
7r aX it aX 

where 

/ n (A) = (l + e*)- 1 . (115) 
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Employing Q114JI and taking into account that e n (— A) = e n (A), the entropy Q91|) can be written as 

'en(O) 



9 /" e n(oo) / 

-V(-l) 5 ™ / de n (/ n (A)ln/ n (A) + (l-/ n (A)ln(l-/ n (A)) 



n=0 



47 -1 

— £(-1)** - L(/M)) , (116) 



7T 

n=0 



where /i, ^ = / n (0), /n 00 '* = / n (oo), and we used the Rogers dilogarithm defined as 

K^-ifd.^+ya). air) 

This function satisfies the following relations 

L(x) + L(l - x) = L(l) , L(x 2 ) = 2L(x)-2L(^—), (118) 

1 + x 

the second of them is called the Abel duplication formula. Moreover, the following identities hold 
(the first is an easy consequence of the Abel formula; for a proof of the second one see, e.g., |16| ) 



y L( _L) + 2 L(-) = L(l), yi( )=-^-L(l), (119) 

k=2 k=l Sm V „+3 I 

where n > 1 and i/ > 2 are integer. 

The knowledge of the entropy allows one to find the corresponding specific heat, C u = T8ts(T). 
On the other hand, at low temperature, we have jSZI, C u = ^-T + o(T), where c is the effective 
central charge of the conformal field theory related to the T = limit, and v s is the speed of sound 
(Fermi velocity), which is v s = \ for our choice of the coupling constant J = 1/tt. Thus, taking into 
account that L(l) = vr 2 /6, we infer from (|116|) that for the alternating chain we have 



k=l 



where we introduced a new variable, f v (\) = 1 — /o(A), and used the first property in l|118|h Observe 
that the replacement of the quasi-particle corresponding to the negative parity by a formal k = v 
particle makes eqs. ()93|) to look uniformly (since e v = — eo, f] v = l/^o)- Furthermore, the formal z^-th 
particle contributes to the central charge (|12()|) in the same way as a genuine quasi-particle. This 
allows us, in the following analysis, to treat all the k = 1, . . . , v particles on the equal footing. 

For T <C 1, the functions ln(l + r] n (X)) vary very slowly for |7rA| <C — InT and for 1 7r A | 3> — InT. 
Therefore, at H = 0, T <C 1, and A — > oo, eqs. (|93[l turn into the following system (which is the same 
as in the homogeneous XXZ case |14j ) 

V 

4 -1 -i= l[(f m r* Ikm , k = i,...,u, (121) 

m=l 

where fk = fi , and ifc m is the incidence (adjacency) matrix of the graph (that coincides with the 
Dynkin diagram of D n with n = v) 
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v — 3 v—2 v—\ 



Under the restriction < / n ,(A) < 1, the system (|121 jl has a unique solution: 

ft^JkTTf' k<u-2, /(-)=/(«) = !. (122) 
It follows from the first relation in (|119|) that 3 



Y,L(fl oo) )=L(l). (123) 

k=l 



A remark is in order here. In the thermodynamic Bethe ansatz approach |3J to the simply laced 
affine Toda models (when the minimal part of the scattering matrix is considered) , the corresponding 
effective central charge c[g] is given by (|120|) . where = and the values of satisfy the system 
(|121|) with Ikm being the incidence matrix of the related Lie algebra q. From this point of view, 
equation (|123f) represents the fact that c[D n ] = 1, independent of n. In the case of Q = A n , the 
system (fT2TJl is solved bv [TUBS] 

2/ 7T 



sin ( — 



h = . , k = l,...,n, (124) 

n+3 



which, according to the second identity in (|118|) . implies 

cK = — — , 125 
n + 3 

as is well known for the A n affine Toda models. 

Now we analyze equations (|93|) for T< 1 and A — > 0. Consider first the case when $i < $2 < v—1. 
The presence of the s(X)5 n ^ i terms in (|93|) leads to e n /T = —00 for n = #i,i?2- Thus, we have 
= fj^ = 1, which implies the following modification of the system (|121|) 

V 

{f£ >) r 1 -i=*k n (/i 0) r |J * m , k=i,...,u, (126) 

m=l 



where a k = (1 - tfjWjX 1 - <5fctf 2/ 

The system (j!26|) corresponds to a graph which is the -D^-type diagram with the $1 and $2 
nodes removed (notice that these nodes contribute cq = 1 each to the central charge). If $2 < ^—2, 
the corresponding graph consists of three disconnected subgraphs - that of A n with n\ = — 1, 
n 2 = — 1) and D n with n = v—'&i. The contribution of the latter subgraph is L(l) and it 

cancels in (|12()j) against (|123|) . Taking into account p25jl . we conclude that the corresponding central 
charge is 

c(v\# u 2 ) = cL^-i] + cL4 w _ x ] + 2c + c[A,_ft,] - c[D v ] = j-^ + , (127) 



3 Which means, in view of 1118H and the relation /„ = 1 — /o, that the pseudo-energies at A = 00 do not contribute 
to the central charge. 
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which agrees with our numerical results (|72[) and with the conjectured formula [20] for the XXX 
alternating chain. For the special case where (Si,^) = (|, 1) we recover the result of j^j. Observe 
that (|127l) remains valid also for <5t? = 0, although the derivation should be slightly modified. 

If ^2 = v— 2, the emerging D n -type subgraph for n = 2 is degenerate (it consists of two discon- 
nected nodes). But we still have c[X>2] = 1 because in this case eqs. (|126|) yield fj^i = fu°^ = \ and 
the contribution of these nodes is L(l) thanks to the fact that L(^) = ^L(l). Thus, eq. (|127|) remains 
correct in this case as well. 

For $2 = eqs. (|93|) imply that €q/T = +oo and hence = ff^ = fjP = 1 (so these 

nodes contribute 3co = 3 to the central charge). Thus, we have to consider the system (|126j) with 
<Tfc = (1 — — 5fc jl/ _i)(l — 5ku)- The corresponding graph consists (for 5$ > 0) of two A n dia- 

grams with n\ = 1 and 712 = v—"Q\—2 = 1. Therefore, we conclude that the resulting central 
charge is given by the same expression (|127(l . Again, if <5t?=0, the derivation should be slightly 
modified but (|127|) remains valid. 



4.4 Remarks on the underlying CFT 

Eq. (|127|) can be interpreted as the sum of the central charges of two SU(2) WZW models at levels $1 
and <5$ which indicates that the structure of the effective conformal field theory is WZW^ 1 (g) WZWs& ■ 
Furthermore, it was argued in for the homogeneous XXZ chain that its underlying CFT contains 
composite fields formed by products of the Gaussian and the Z(t?)-parafermionic fields in agreement 
with the representation c(z/|i?) = 1 + c[^_i] (cf. (|127|) for <5$ = 0). We may expect a similar 
operator content for the alternating chain, in particular, composite fields involving Z{$i) and Z{5 r &) 
parafermions are likely to be present. 

It should be remarked here that (J127|) gives the resulting central charge but does not reflect the fact 
that all v quasi-particles in (|92|) contribute to its value. It is however known for the thermodynamic 
Bethe ansatz treatment of models with Lie algebraic symmetries, like the simply laced affine Toda 
models, that the related Lie algebraic structure is often exhibited in the quasi-particle content of 
characters of the related CFT [431 144) . It is therefore natural to look for a basis in the CFT related 
to the (alternating) XXZ chain that would contain v quasi-particles. 

It was observed in j^SJ that many characters of various conformal models admit the following 
(so-called fermionic) representation, 



x(q) = /~2 Q 



ira'Bm+m-a+no 



mi , 



n 

fc=i 



((l-S)m + u) fc 



(128) 



where B is nxn real matrix, a and u are n-component vectors, ao is a constant, and the (/-binomial 
coefficients are defined for < m < n as 



(g)n 

{q)m{q)n-r, 
1 

{q)m 



if 
if 



n < 00 
n = 00 



(9), 



IK 1 

k=l 



(129) 



Usually B is related to the Cartan matrix C g of some Lie algebra q, and there are typically certain 
restrictions on the summation over m k . Analysis of the q — > 1~ asymptotics of ()12<S[) allows one 
to find the corresponding central charge as a combination of Rogers dilogarithms with arguments 
satisfying two systems of non- linear equations 43 a . It turns out that our equations ()120j) , (|121|) , and 
()126j) are exactly of this form provided that we make the following identification 



n 



B 



00 . 



U u -1 



(130) 
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where Cp v is the Cartan matrix of the Lie algebra D u . Furthermore, characters of the type ()128(1 
based on q = arise in the context of the spinon basis of SU(2) WZW models related to the 
homogeneous XXX spin chain 45 . This, together with the correspondence (|130f) . makes it plausible 
that the spinon basis related to the (alternating) XXZ spin chain with the anisotropy 7 = ir/v will 
correspond to ()128|) based on 5 = D u . 

Let us stress here that the correspondence (I13U() is based solely on the data obtained by the 
thermodynamic analysis of the ground state of the system. It is worth to mention that an additional 
information, in particular, about a and ao in ()128j) . can be obtained from the finite size analysis 
of low lying excitations (see, e.g., |53[ I17j) because this analysis allows one to find the anomalous 
dimensions of the fields present in the underlying conformal theory. 



4.5 The H = 0, T » 1 limit 

It is also interesting to examine the high temperature behavior of the entropy and deduce the number 
of states of the model. In the high temperature limit the pseudo energies e n become independent 
of A Jl] . Consequently, neglecting the free terms in (|93|) , we obtain the same constant TB A equations 
as in (jnij>. Thus, we have (cf. (JT22J)) 

1 + Vn = (n + l) 2 -^-i , n = l,...,u-l. (131) 

According to (|94|) and (|95|). the free energy density at H = in the high temperature limit becomes 

f = ~J Yl (1 + <W-i) m(l + . (132) 
Hence, taking into account (|131j) . the entropy is given by 

5= Y E Hn + l). (133) 

Which implies that (recall that i?j = 2Si ) the number of states for the system is 

((2Si+l)(2S 2 + l)) T , (134) 

as one could have expected. 



5 Generalization to I spins 

Consider the XXZ chain containing I different spins, Si, ... ,5*/, assigned to sites with periodicity I, 
i.e, the spins at the sites n and (n+l) are the same. Let us outline how our analysis of the alternating 
chain extends to this case. 

Introduce the following set of transfer matrices for i = 1, . . . , I 

r» (A) = tr a (it$$ (A) B**zl (A) . . . R S a f + \ W <f (*)•■■ W ^ff (A)) . (135) 

They commute with each other and, extending the proof given in Section 2, it can be shown that 
n!=i r ^(0) generates a shift by / lattice sites (see Appendix B). This means that the Hamiltonian 
constructed as follows 

T 1 

n = q^d x lnr^(\) a=q (136) 

i=l 
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defines (in the gapless regime) at zero temperature a conformally invariant system. We remark that 
Ti is local but involves now interactions of (1+1) nearest sites. 

The Bethe ansatz equations now read 

— / sinh7(A Q + iSj)\ N ^ 1 -^4 sinli7(A Q — Xp + i) 



n(^tm) -n s^ . -i.-.^ <-) 



Denote i?j = 2Si. For the Bethe ansatz analysis, the order of the spins on the lattice is irrelevant, so 
we assume in the following that $i < i?2 < ■ ■ • < $i < v. Then, at T = 0, H = 0, the ground state in 
the thermodynamic limit is given by Pfc(A) = js(\)5k$ i (the Dirac sea of Z different strings) with the 
corresponding energy 



l l 
j7r Yl /^^(A)KA) = ^^[^(^ + i)-^(^ T ^ + i)]. (138) 

A — 1 J A A — 1 



60 = -jr 



i,i=l_oo »,i=i 
At T > the free energy at equilibrium is given by 
l 

i » 

/ = e 



tE{/ dAs(A)ln(l + %i (A)) + ( ^_ 1 / dA S (A)ln(l + 7 ? - 1 (A))}. (139) 

* ^ •/ — oo J — oo 

Equations (|93|) for the pseudo-energies are modified only in their free terms by the replacement 
°f I Si=i w ith j Yl\=i ^ktii • Observe that this change does not affect the A — > oo analysis in the 
Section 4.3. For A — > we have now f$ i = 1 for i = 1, . . ., I. That is, removing Z nodes, we obtain A n 
diagrams with n\ = — 1, ri2 = $2 — — 1, etc. This yields the following central charge 

cH^i, ...,*,) = + £ , Mi = tf m - tf, . (140) 

Notice that the condition i?j < ^ implies that Z < i/— 1. It is easy to show that, for a given v, the 
maximum of (|14U|) is attained if Z = i/— 1 (that is, if all the possible spins are present), 

maxc(i/|{tfi}) = c(i/|l,2,... 1) = v-l. (141) 

On the other hand, if Z is fixed, then the configuration i?j = i corresponds to the minimal possible 
value, c = /. 

The analysis carried out in Section 4.2 extends to the case of Z spins by considering the /-component 
counterpart of the Wiener-Hopf system (|1U5|) (with £ = u /(v — di) and the coefficients in front of e _7rA 
being hjl). Then analogues of Q1U8|) and (jlU9f) hold (with a factor 2/Z appearing on the r.h.s. of I|1U9|) ). 
and we derive tfS^m) = ±(G+(0)G-(p)) u . Observing that (G+(0)G_(0)) H = i(0)^, we obtain 
k*S(m) = . Using now ()139|) and proceeding as in (|lllj) . we conclude that 



Whence we infer that 



/ = ,_iJ*a». (142) 

27T 1/ — 



-x-f-^V ("3) 



that is, the magnetic susceptibility x depends, for a given anisotropy, only on the maximal spin Si. 
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6 Conclusion 



In this article the alternating quantum spin chain model with periodic boundary conditions and 
its Z-spin generalization were investigated via the finite size effects and the algebraic Bethe ansatz 
method. In a previous work |26| the boundary free energy for the alternating chain was computed 
by means of the thermodynamic Bethe ansatz method. On the other hand, the effective central 
charge of the open XXZ spin chain with general diagonal and non-diagonal boundaries was computed 
recently by means of finite size corrections |48| . A similar computation, but for diagonal boundaries 
only, was also undertaken in |49j . In both cases the effective central charge clearly depends on the 
boundary parameters of the model. Let us also mention that the finite size correction approach, when 
boundaries are present, provides somehow the same information with the TBA description in the so 
called R channel [HOI , whereas the thermodynamic Bethe ansatz approach see, e.g., [13 021 EH] , can 
be compared with the TBA description in the L channel |491 152j . We hope to generalize the finite 
size correction approach for fused models with open boundaries. 

To complete the picture, the thermodynamics of the alternating spin chain and the RSOS^i, $2; q) 
model with integrable boundaries should be also investigated. Finally, it would be of great interest to 
continue the thermodynamic investigation for spin chains and RSOS models related to higher rank 
algebras. We hope to address these questions in a future work. 
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Appendix A: universal XXZ R-matrix 

The universal R-matrix (JHj was given (in a slightly different form) in without an explicit deriva- 
tion. We find it instructive to provide here a proof of l[H|). Technically, it extends the derivation used 
in the case Si = S% [HI (see also [2E1 EH EH ) • We will also prove some useful related formulae. 
Recall that the co- multiplication on U q {sl2) is defined as a linear homomorphism such that 

A(S ± ) = S ± <g> q~ s3 + q s3 <g> S ± , A{S 3 ) = S 3 ® 1 + 1 <8 S 3 . (144) 
Evaluating (JBJ) in the representation V\ ® Vs 1 <8> Vg 2 (and interchanging A with fj,), we obtain 

Rf^iX) L S 2 \X + fj,) Lf(/i) = Lf (/x) Lf 2 (A + M ) i?f 2 lS2 (A) (145) 

which imposes the following conditions on R lS2 (X): 

R{X)A' = AR(X) , (146) 
R(X) (e 2e7A 5 e <g) q~ s3 + q 3 * tg> S £ ) = {e 2e ^ x S e <g> q s3 + q~ s3 <g> S e ) R(X) , (147) 

here (and below) e = ±, and A' is the opposite co- multiplication, i.e., = A l / q . These equations 
determine the universal R-matrix uniquely (up to a scalar factor) [HI EE]- Denote R± = (i?(±oo)) _1 . 
Then the asymptotics of (fli5 |) -(fl47 |) yields 

R e A = A'R e , (148) 
R e (S e ®q S3 ) = {S e q- s3 )R e , R e (q~ S ' 3 ® S~ e ) = (/' ® S~ e ) R e . (149) 
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These equations imply that R± are the constant universal R-matrices (jlUj) . Observe that (1148(1 
implies that R e maps the A-basis in Vs 1 ® Vs 2 into the A'-basis. That is, for the orthonormal set of 
eigenvectors \j,m), m = j,j — 1, — j of a projector Vj we have 

R^\j,m)=^\j,my , (150) 

where the scalar factor <p e - does not depend on m (its explicit form is given below). 
Now we introduce r±(A) = R(\)R± and then (|ll£j) - ((TiT|) acquire the form 

r e (A)A = Ar £ (A) , (151) 
r e (A)(e 2e7A 5 e ®g s3 + ( J R e )- 1 (g s3 (g)5 e )i? £ ) = (e 2e ~< x S e ® q s3 + q~ s3 ® S e ) r £ (\) . (152) 

It follows from (|151f) that r =t (A) are functions of the spin operator J, i.e., 

5i+5 2 

^(A)= £ rf(X)V 3 , (153) 
|Si-S a | 

where ^(A) are scalar functions and Vj is the projector onto Vj in Vs 1 <8> Vsr 2 . 

Using (|144|) . (|149|) . and (|150|) . it is straightforward to check the following relations 

[S € <g> q s3 ,A{S e )} = [q- s3 ® S e , A{S e )} = [(RerHq S3 ® S e )i? e) A(S € )} = , (154) 
(S* ® 9 s3 + 9< 2+2 ^9- 53 » S e ) |j, ej) = (10 9 253 ) A(5 e ) |j, ej) = , (155) 
(V 3 S*)R e + q< 2+2 ^R e (S* ® / 3 )) |j, ej) = (q 28 * ® 1) A'(S e ) |j, ej}' = . (156) 

Eqs. (HED-JEU) show that we can apply (|152|) to a highest/lowest weight vector |j, ej). Doing so 
and taking (|155|) - (|156|) into account, we obtain 

rj+iW (e e7A - q < 2+2 ^ e -^ x ) = (e e ^ x - q -<^i) e ~^) r *(\) . (157) 

Hence 

^ +1 (A) = -9-^)|^|^(A). (158) 

Solving this functional equation and fixing the overall normalization, we obtain ©. 

Notice that Q has the form g^ e ®f2(A, 3)(R e )~ 1 , where the functions i/j e ( x ) an d f2(A,x) are in- 
variant with respect to the replacement of q by q' 1 , and, moreover, U(X, x)Q(— A, x) = 1, ip- e (x) = 
—i/j e (x). Whence, taking into account the property (R e )i/ q = (i? e ) _1 , we derive the first symmetry 
relation in (|13|) 

R 1/q (-\)=q- Mr) n(-\J')R e =R t q-^(n(\J)y 1 = (^(A))" 1 . (159) 

Here J' = 3x/ q is the two-node spin operator corresponding to the opposite co- multiplication and we 
used that R e J = fR e , as follows from (|148|) . 

If the algebra generators are represented in terms of matrices as in (J3J, then it is easy to check 
that J* = J and (R+Y = (R-) -1 , where t stands for the matrix transposition. Whence we derive the 
second symmetry relation in (|13|) . 

(R(-X)f = R_ 6 q^n(-XJ) =R_ t q-^W(n(\J)y 1 = (^(A))" 1 . (160) 
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It is straightforward to see from (|15|) - (|16|) that action of the *-operation (|38j) on the two-node 



spin operator 1 is given by (JT g ) 
the first relation in (|3*9*)): 



Jg. This, along with the property {R±, q 



(R 



(R{\)Y = R^ jq n(-XJ q )q^-^ = (R q -(X)Y 



leads to 



(161) 



The second relation in (|39j) follows from (|22|) . (|161j) . and the symmetries 1)1 3jl and (|23fl . For instance, 
for 7 E R we have 



-1 „7A(1®S 3 ) 



-7A(1®S 3 ) ^(_A) e 7A(l®53) = j^(_~ x) = (^(A))- 1 . 



(162) 



To prove the important identity (jllj) . we employ Theorem 3.1 from |47| that states (in our nota- 
tions) that the U q {sl2) Clebsch-Gordan coefficients satisfy the following relations 



r £1 ^2 j 

mi tti-2 m 
mi ttt-2 m 



(•_ 1 \Si+5 2 -i g «(Si)+«(5 2 )-«(i) 
(_l)-Si+S 2 -i 



S 2 Si j 

777,2 Wl 777 



£2 Si j 

7772 mi TU 



(163) 
(164) 



where k(x) =x(x + 1), and (^+ lS2 )ml'm2 = ( m i\ ® { m 2\R+ lS2 \ m i) ® \ m 2) ar e matrix elements of the 
R-matrix in the Vs 1 <S> Vs 2 basis. These formulae show that for (|15U|) we have 

„«(Si)+«(S 2 )-«0') 



and also that 



(mi) (m 2 | J R+ ,S |j,m) = (-1) 2S ® ("4 1 F |j> m } , 



(165) 



(166) 



which proves (|lljl for e = +. Analogues formulae for e 



are obtained by observing that the 



map S H 



S i $ 



-S 3 , q «-> g 1 is an automorphism preserving the co- multiplication (|144|) . In 



particular, we have 



which along with (|15Ujl and (|165j) leads to the following relations 



Inverting these relations, we obtain 

Vr 52 (^) 2 -^e 



S1+S2 

R+R- = £ C^) 2 ^ 

j=5S 



(167) 



(168) 



n 



_e(«(Sfx)+K(S 2 )-«0')) 



(169) 



Indeed, to prove this formula, it suffices first to observe that, as seen from (|168|) . the only projector 
which is present in all the factors of the r.h.s. of ()169j) is Vj, and then to check the corresponding 
normalization. In the simplest case, (Si,S2) = {^,S), eq. (|169j) yields 



is 



'[25 + 1], 



fC*S+i)_ q ( R hSyi R kj 



(170) 
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Consider now the universal R-matrix introduced in (|22jl. It is a solution to equation (|145|) for 
the L-operator 1)20(1 . The latter, unlike the L-operator (|4"]l. does not decompose into two Borel 
components. As a consequence, instead of (|146j) - (|147j) . we have four A-dependent relations: 

R(X) (e^ /X S e ® q~ sS + q s3 <8> 5 e ) = (e £7A S e ® + g" 5 " 5 e ) i?(A) , (171) 
R(X) (e^ x q- s3 <g) S e + S e ® g 5 ') = (e £7 V 3 ® S e + S £ g> ^ 3 ) A(A) . (172) 

Observe that R(X) = (R(X)f solves the same set of equations, which, taking into account that the 
solution to (|171[) - (|172|) is unique (up to a scalar factor), implies that (R(X)) 1 = ((X)R(X), with ("(A) 
being a scalar function. In fact, £(A) = ±1 because tot = id. The possibility of £(A) = —1 is excluded 
by the fact that the vector \S±) ® IS2) is an eigenvector of R Sl 2 (A) with a non-vanishing eigenvalue 
(hence the first matrix element of R SlS2 (X) in the basis corresponding to (J2J) is non- vanishing) . Thus, 
we infer that C(A) = 1, which proves the last relation in (|2Mj) . The remaining symmetries in ({SB} are 
easily derived from this symmetry and the above established relations (|13|) . for instance 

R(—X) = e -rA(i®^ )i2( _ A)e7 A(i®s3) (173) 

= ( e -7A(l®53) (i?(A))t e7 A(l®S3))-l = ((^(A))*)- 1 = (^(A))- 1 . 

Finally, since R(0) = R(0), the symmetries l)14|l follow obviously from IJ13JI and (|23|) . 



Appendix B: Construction of the shift operator 

Let Z and iV be integer such that I divides N. Introduce the following invertible operators 

N N 

q i i= n <ft-^(°)' &*= n R i%-m^ ( i74 ) 

fc=i mod I k=i mod £ 

S- S ■ 

where Rn,m (A) is the XXZ R-matrix (assigned to the sites n and m of the lattice in the sense of 
Section 2). Notice that in each product all the factors commute with each other iff i ^ j mod I. 
Along with the property (|45|) this implies that 

Q% j Qt* = Q 3 t = 1 0JV • (175) 
Consider now the transfer matrices (|L-S5[) . It is not difficult to derive for them an analogue of eq. I|27jl. 

T (i) (0) = u(i) qM-1 qM-2 _ _ _ gi.1 Qi,l Qi,l-1 _ _ _ ^ (176) 

where Z/W = Yl k=i mod i ^fe,fc+£ (the- factors are ordered from the left to the right) acts nontrivially 
only on the nodes n = i mod I, 

£ n = £ n+/ if n = imodl. (177) 

Observe that U = U^U^ . . . is the lattice shift operator (generating a shift by I sites). 
Notice that Q^U^ = U^Q^. Therefore (taking all to the left) we have 

rW(0)r( 2 )(0)...r«(0) = U Q 1 ! 1 Q 1 ! 1 ' 1 . . . Q 1 / Q 1 ? (178) 

x Q^Q 2 - Z ...Q 2 _ 4 Q 2 - 3 

x Q l + 1 > l - 2 Q l + l ' l ~ z ...Q l + l ' l Q l : 1 > 1 
x Q 1 {- 1 Q 1 J- 2 ...q¥q 1 ± 1 . 
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Notice that here we have Q_ if i < j and if i > j. Let us show that all Q-'s can be led 
to meet their Q + counterparts and thus be canceled by the relation Q175|) . As the first step, the 
rightmost Q_ in each row, i.e, Q l ^ +1 , i = 1, . . . , l—l, is canceled against the leftmost Q + in the next 
row. Observe that this eliminates obstacles for moving then the second from the right Q_ through 
the entire next row (obviously, Q l ± commutes with Q± J if i,j,i',f are all distinct). Thus, at the 
second step, Q l ^ +2 for i = 1, . . . , 1—2 are canceled against Q % + 2 ' % (which, after the first step, became 
the leftmost factors in their rows). Proceeding this way, we eliminate at the fc-th step l—k pairs 
(Q l - +k , Q l + k ' 1 ) and also make possible the next step (that is, Q^ +k+1 w [\\ De to move through 
the next k rows). Consequently, after the (I— l)-th step we will have all Q's eliminated, and so we 
conclude that t^(0) . . . r^(0) = U. 
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